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ABSTRACT 

A quasi-steady state model of the consequences of rotation on the hydrodynamical 
structure of a stellar radiative zone is derived, by studying in particular the role of 
centrifugal and baroclinic driving of meridional motions in angular-momentum trans- 
port. This nonlinear problem is solved numerically assuming axisymmetry of the sys- 
tem, and within some limits, it is shown that there exist simple analytical solutions. 
The limit of slow rotation recovers Eddington-Sweet theory, whereas it is shown that 
in the limit of rapid rotation, the system settles into a geostrophic equilibrium. The 
behaviour of the system is found to be controlled by one parameter only, linked to the 
Prantl number, the stratification and the rotation rate of the star. 
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1 INTRODUCTION 



The study of the effects of rotation on the hydrostatic struc- 
ture of a star is a long-standing problem (Eddington, 1925, 
Sweet, 1950). It is the key towards a better understanding 
of stellar evolution, and consequently, of the whole observ- 
able universe. The effects of rotation on the stellar evolution 
can be separated into two classes: effect on the contraction 
or expansion with time of various stellar regions and rota- 
tional mixing (of chemical species and angular momentum) 
through waves, instabilities and centrifugal driving of merid- 
ional motions. The effects of rotational mixing in stars was 
reviewed by Pinsonneault (1997) and can be directly ob- 
served (e.g. Li abundances and Li dip problem, MS turnoff, 
low-mass giants deep mixing ...). Recent works by Maeder 
& Meynet (2000), or Chaboyer, Demarque & Pinsonneault 
(1995) for example, attempt to consider the effects of ro- 
tational mixing on stellar evolution through the resolution 
of one-dimensional stellar models using various ad-hoc pa- 
rameterizations. These models are successful in explaining 
some aspects of the observations but their predictive power 
is sometimes limited. 

In this paper I propose to approach the problem through 
a new path, which has different limitations and advantages 
thereby providing an interesting complement to the existing 
models. I choose to study the nonlinear effects of rotationally 
driven meridional flows on the angular-momentum distribu- 
tion of a stellar radiative zone. In the limit of a slowly ro- 
tating star, this is equivalent to the well-known Eddington- 
Sweet problem. In the case of rapid rotation, this problem 
has received less attention; it can now be studied through 
the formalism introduced in this paper. This type of ap- 



proach has been considered in the past (e.g. Tassoul & Tas- 
soul (1983)); however, in previous works the nonlinearity 
of the angular-momentum transport processes was usually 
neglected mainly by assuming that the star was rotating 
nearly uniformly. As I shall show, this assumption is not 
self-consistent. 

In Section 2, I present the equations for the hydrody- 
namic structure of a quasi-steady, laminar and axisymmetric 
radiative zone. The numerical results suggest a new scaling 
of the variables of the problem which can then be solved, 
in certain limits, analytically. The solutions are presented in 
Section 3. They provide the first self-consistent model of the 
effects of rotation on the centrifugal driving of meridional 
motions (assuming that the flow is a quasi-steady, laminar 
flow) and relies on one parameter only: the mean rotation 
rate of the star. The results are discussed in Section 4. 



2 THE MODEL 

I shall consider solar-type stars only and study their quasi- 
equilibrium structure under rotation. The overlying convec- 
tion zone is affected by rotation mainly through the Corio- 
lis distortion of the convective eddies. These effects are ex- 
tremely complex; they are simply represented in this work 
through the rotation profile imposed by the convective zone 
to the underlying radiative zone. The radiative zone is as- 
sumed to be a stable, isotropic fluid with uniform^] dynam- 
ical viscosity jj,. 



* This is a fairly good approximation; in any case, the viscosity 
only affects the dynamics of the system in thin boundary layers. 
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The steady-state structure of the rotating radiative zone 
is obtained by perturbing the hydrostatic equilibrium in the 
following way: 

p h u ■ Vu = -Vp - /5 h V$ - |5V$ h + mV 2 u + ipV(V-u) , 

p h T h u ■ Vs h = V-(fcVT) ,£. = _£. + _ 
Ph ph Jh 

V 2 <E> = 47rGp , V-(phu) = . 



(1) 



where p is the pressure, T is the temperature, p is the den- 
sity, $ is the gravitational potential, u — (u r ,ug,u^) is the 
velocity field, s is the entropy and k is the thermal con- 
ductivity. As the system is axisymmetric and near-uniform 
rotation is not assumed, there is no need to work in a rotat- 
ing frame. Quantities denoted with suffix h are derived from 
the non-rotating hydrostatic equilibrium equations, and per- 
turbations to that state are denoted with tildes. The equa- 
tions are linearized with respect to the background state 
(e.g. 4> <g $h, p <C ph), but the full nonlinearity is kept in 
the flow and in particular in the advection process u ■ Vu. 
This approach is rigorously correct provided the star is far 
from breakup (that is, when the mean ellipticity due to rota- 
tion is much smaller than unity). The steady-state assump- 
tion is valid provided the dynamical timescale of this system 
is short compared to the stellar evolution and spin-down 
timescales. This is not always necessarily the case, so that 
the direct applicability of the results presented in Section 3 
must be checked a posteriori. 

These equations must be supplemented with appropri- 
ate boundary conditions. Strictly speaking, correct bound- 
ary conditions would require a theory of the rotation and 
circulation in the convection zone to which the interior flows 
must be matched. Here, I shall simply assume that the 
forces maintaining the flow in the convection zone are so 
high that the angular velocity is unaffected by the flow in 
the underlying radiative zone, and is given by f2 cz (r, 9) = 
X^^Lo ^™ cos ™$ (where 9 is the co- latitude and r is the radial 
variable). In addition, I suppose that the stresses affecting 
the meridional circulation are viscous-like. Since the effec- 
tive viscosity in the convection zone is much greater than 
the viscosity in the radiative zone, it is not implausible that 
the continuity of stresses across the boundary imposes the 
constraints Q, — J7 CZ , r r ,g = rd(ue/r)/dr + (1 / r)du r / 89 = 
and 80,/dr is continuous (and therefore null) at the inter- 
face. The inner core (for r < ri n ) is removed from the com- 
putational domain in order to avoid singularities. This core 
is impermeable and is assumed to rotate rigidly with angu- 
lar velocity f2j n determined in such a way as to ensure a null 
angular-momentum flux through the boundary. Finally, the 
regions outside the domain of computation are assumed to 
satisfy V 2 T = where T vanishes both at r = and when 
r — > oo (since the inner core supports no fluid motion and 
the effective thermal conductivity in the convection zone is 
assumed to be much larger than in the radiative zone). 

Using the assumption of axisymmetry, I reduce the mo- 
mentum equation in (1) to: 



u ' V{ (£ sin9u^) — E U T> (£ sm9u^) , 

1 d / 2 \ _ sin6> / 8p h <9<£> 1 dp 

+ E V D 2 (£ sin0av) , 



where l> = Vxu is the vorticity, £ = r /r c is the new normal- 
ized radial coordinate, z is the normalized cylindrical coor- 
dinate that runs along the rotation axis, E v — v/r^Sl* is the 
Ekman number (where v = p/ ph), and e = r- 2 SlJ(9$h/90 _1 
is the ratio of the centrifugal to gravitational forces. The first 
of these equations simply describes angular momentum con- 
servation, relating the advection term (LHS) to the diffusion 

term (RHS)pJ the second equation is the azimuthal compo- 
nent of the curl of the momentum equation, commonly re- 
ferred to in the geophysical literature as the 'thermal wind 
equation'. In all these expressions the following normaliza- 
tions have been applied: [r] — r c , [u] — r c Q,+ , [$] = r^fi 2 ., 
[T] — IK, [p] = lg/cm 3 , where r c is the radius of the radia- 
tive zone and Q* is the typical rotation rate of the star. The 
operator D 2 is defined as 



D 2 = 9 



sm9 8 

+ 



V sin6» 89 J 



(3) 



The energy equation becomes, to first order in the ther- 
modynamical perturbations 



0,1 E v 



(4) 



where a — pc p /k is the Prandtl number, and is the 
background buoyancy frequency. 

Finally, the equation of state can be combined with the 
radial and latitudinal components of the momentum equa- 
tion to provide an expression for p: 



±21 = E^rlnl 

ph 89 p h 



cos# 

1 

shit* 



<9$ 



%w (5) 



This expression is used in the thermal wind relation in (g) as 
well as in the latitudinal derivative of the Poisson equation 



|-V 2 $ = 47rG 
89 



Ph 



cos# 2 



9$ 



PhdT 
T h 89 



(6) 



The resulting equations (^|), (^) and ([]) combined with 
the mass continuity equation are solved for the unknowns 
u r , ue , u,f„ 9$ j89 and T. 



3 NUMERICAL AND ANALYTICAL RESULTS 

The numerical study of the system presented in Section 
^| was performed using the standard solar model calcu- 
lated by Christensen-Dalsgaard et al. (1991) as the refer- 
ence state. The imposed latitudinal shear was chosen to 
be that observed in the solar convection zone (with Q cz — 
1 — 0.15 cos 2 — 0.15 cos 4 # in units of O*). As the rotation 
rate fi* is varied, it appears that only two parameters con- 
trol the behaviour of the system: the Ekman number (which 
controls the boundary-layer behaviour and the flow veloci- 
ties) and a new number 



A = o-iVh/fi 2 



(7) 



(2) 



which controls the dynamics of the bulk of the radiative 
interior. 

T The diffusion term is expressed here in the case of isotropic 
viscosity, but it is worth noting that any other Reynolds-stress 
prescription for anisotropic transport can also be used. 
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3.1 Slow rotation case: A ^> 1 

In the case of slow rotation, the numerical results suggest 
that T and the poloidal components of the velocity u r< e scale 
with E v and A as: 



T = eT^T , u r ,e = E v j\u r fi 



(8) 



where the quantities with bars are the scaled quantities, of 
order of unity. It is also found that and $ are always of 
order of unity, which is expected. Note that the scaling for 
the meridional motions is a local Eddington-Sweet scaling 
(see Spiegel & Zahn, 1992) . Using this ansatz into the system 
of equations given in (1), an expansion in powers of 1/A 
reveals that the angular-momentum balance is dominated 
to zeroth order by viscous transport through: 



D 2 (£ sme U4> ) = 



(9) 



which determines the angular velocity profile uniquely. Us- 
ing this result in the first order equations provides a relation 
between the temperature and gravitational potential pertur- 
bations: 



T d ( 1 

az \T h 



^-v 2 4> 
de 1 



sine | ^ 



4ivGp h r c 
9h 



dlnp h 



d£ d6 

cose 9 



(10) 



sin# 



9$ 

~de 



dT 



which can be solved independently for T and d^/de. Fi- 
nally, the temperature fluctuations lead to meridional mo- 
tions through the energy advection-diffusion equation: 

iV = V|T. (11) 

In Fig. I show the results of the numerical solutions 
for the angular velocity profile and the meridional motions 
corresponding to a slowly rotating solar-type star (for which 
A ~ 10 4 ). The angular velocity profile is exactly the same as 
one obtained through solving (V 2 u)^ = 0. The meridional 
flow velocities decrease strongly with depth and the stream- 
lines follow a layered structure reminiscent of a Holton flow. 




Figure 1. Numerical solution for a solar-type star rotating 100 
times slower than the sun (A ~ 10 4 ). The quadrants show the 
radiative zone only. The upper panel shows the angular velocity, 
which is viscously dominated. The interior rotation rate is 0.957 
times that imposed at the surface at the equator. The lower panel 
shows the streamlines (dotted lines represent a clockwise flow, 
and solid lines represent an anti-clockwise flow). The contours 
are logarithmically spaced. 



3.2 Rapid rotation case: A<1 

In the case of rapid rotation the correct scaling seems to be 



T = XeThT , ibr,e = E v u r . 



(12) 



This time, I perform an asymptotic expansion in the small 
parameter A. In this limit the temperature fluctuations are 
strongly damped by the rapid heat diffusion (as A <C 1 is 
equivalent to the small Prandlt number limit) and the sys- 
tem reaches an equilibrium that is determined by the zeroth 
order equations: 



sin# 



d£ He 



mse 2 



(13) 



These equations describe a geostrophic equilibrium, and can 
in principle be solved for u 2 ^ and $ alone. The solutions pro- 
vide, to the next order in A, an equation for the meridional 
flow through the advection diffusion balance: 



u ■ (£ sin#u^,) = D (£ sixiOu^) 



(14) 



and, finally, the temperature fluctuations through 



(15) 



The results of the numerical simulations for small A (A ~ 
10 -2 ) are shown in Fig. |^. It is interesting to note that 
the radial shear is much larger than the latitudinal shear, 
which suggests to consider dQ/d£, ^> dfl/de to perform a 
first analytical approximation. Similarly, I suppose that if 
* = d$/de then dty/d£ > 9* /de. These approximations 
can be used into equations (jij) and ([]) to yield the system 



^ di y di 



df\ 4nGr 2 pl, r r „ 
_ cf2L rq, + K \ 

Ph 



(16) 



where K — cos#£ 2 f2 2 , which can easily be solved numeri- 
cally for K and However, if the radiative zone is crudely 
approximated by a polytrope of index 3 then Tk oc l/£ and 
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Figure 2. Numerical solution of the system (1) for a solar-type 
star rotating 10 times faster than the sun (A ~ 10~ 2 ). The upper 
panel shows the angular velocity, which varies strongly with depth 
and the lower panel shows the temperature fluctuations. Note 
that even when the stellar oblatcncss is of order of 10 — 3 , the 
temperature fluctuations remain of order of 10 — 6 . 



Ph oc l/£ 3 between £ = 0.3 and £ = 1. This leads to 



a dlnT h 
-pr and 
£ 



(17) 



Ph ^ <9£ £ ' 

These expressions can be used in the system ( |l6| ) to pro- 
vide a homogeneous relation between and K, which can 
(amongst other power law/oscillatory solutions) be solved 
analytically with ~ —K ~ constant; this corresponds to 
fi(£) oc l/£. The comparison of this very simple analytical 
prediction with the results of the fully nonlinear numerical 
simulation is presented in Fig. |j| in the case where the con- 
vection zone shear is null and in the case where it is of order 
of 30% of the mean rotation rate. When the imposed shear 
is small, the analytical prediction provides a good approxi- 
mation to the numerical results. 



0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 
f 

Figure 3. Comparison of the analytical prediction f2(£) oc l/£ 
(dotted line) with the numerical results (solid line) in the case of 
a solar shear (where the poles are rotating 30% slower than the 
equator) and in the case of no shear for a star with A ~ 10~ 2 . The 
no-shear case satisfies the spherical approximation intrinsically 
well, so that the analytical approximation is closer to the true 
solution in that case than in the sheared case. The discrepancy 
with the true solution increases near the core where the polytropic 
approximation to the equation of state is no longer valid. 



linked with the lower boundary conditions and may not rep- 
resent what can be expected of a stellar radiative zone. In 
particular, the numerical simulations reveal the presence of 
an internal shear layer (related to the internal Stewardson 
layer in the Proudman spheres problem (Stewardson, 1966)), 
which is related to the presence of the impermeable lower 
boundary (and therefore mostly irrelevant to stellar hydro- 
dynamics). This shear layer is reflected also in the temper- 
ature distribution (see Fig. ^). However, as the meridional 
flow and temperature distribution affect the angular velocity 
profile only to the next order in A, it is possible to calculate 
the flow pattern as if there was no lower boundary using 
the simple ansatz Q oc l/£ into equation ([lif). Using the 
polytropic approximation again one can derive 



ip = £ sin 9 



hi 



1 



cos# 



I - < <>s6») 



(18) 



where ip is the scaled stream-function defined as u = 
— Vx(^>/£ sin#). This result suggests that in stars the di- 
mensional radial component of the meridional flow varies as 
U r oc v/r independently of the stellar rotation rate when 
A < 1. Note that this scaling is compatible with the results 
of the work of Talon et al. (1997). As Busse (1981) claimed 
that rapidly rotating systems settle down to purely zonal 
flows in the inviscid case, it is not surprising to find that the 
typical timescale for steady circulation in this rapid-rotation 
limit is the very slow viscous timescale. 



3.3 Lower boundary problem 

The numerical solutions for the angular velocity and the 
gravitational-field perturbation know little of the presence 
of a lower boundary; the meridional motions and the tem- 
perature fluctuations on the other hand seem to be strongly 



4 DISCUSSION AND CONCLUSIONS 

In this paper I analysed the consequences of rotation on 
the driving of meridional motions in stellar radiative zones. 
For a slowly rotating star, this problem has been studied by 
Eddington (1925) and Sweet (1950) who established that 
meridional motions are driven by stellar baroclinicity re- 
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lated to the rotation. For rapidly rotating stars, the problem 
has never been solved self-consistently. I have constructed a 
quasi-steady model of the rotating radiative zone of a solar- 
type (main-sequence) star which can be solved numerically 
for a wide range of rotation rates; through this numerical 
analysis I have been able to show that the effect of rotation 
on the stellar structure can be described with two param- 
eters only: A = aN^/Q^ and E v — i/yVcfi*. The striking 
reduction of this complicated nonlinear problem to a rather 
simple form stems from the very close thermal advection- 
diffusion balance; this is intrinsically related to the very 
small Prandtl number typical of stellar radiative zones. In 
comparison, studies of the effect of rotation on stratified 
shear flows in the high-Prandtl-number limit reveal a very 
different behaviour (e.g. Clark et al., 1971). 

The system approaches asymptotic limits in the case 
of fast rotation, A < 1, and slow rotation, A 3> 1. In both 
cases, interesting analytical reductions of the problem can be 
used. In the limit of slow rotation it was shown that typical 
flow velocities scale as the local Eddington-Sweet velocity, 
and the typical adjustment time of the system to the steady 
state is an Eddington-Sweet timescale (as expected). As a 
result, for molecular viscosities typical of solar- type stars it 
is unlikely that the system would ever relax to the steady 
state found here. 

In the limit of rapid rotation (A <C 1), it was shown that 
the system settles into a geostrophic equilibrium; this can 
be achieved on a rather short timescale, which is linked with 
the generation and dissipation of inertial waves. As a result, 
this type of equilibrium is likely to be present in stars^|. The 
resulting rotation profile varies with radius roughly as 1/r; 
the precision of this approximation depends on the degree of 
shear imposed by the convection zone. Assuming this simple 
profile for the angular velocity, the meridional motions can 
be deduced from angular-momentum conservation. The typ- 
ical flow velocities in the steady-state are of order of u/r c , 
and are too slow to have any effect on material mixing. How- 
ever, the continuous adjustment of the radiative zone to the 
stellar spin-down may lead to some significant mixing near 
the interface with the convection zone through the genera- 
tion and dissipation of inertial waves. This can be studied 
self-consistently only through a time-dependent analysis of 
the flow, which is the subject of future work. 

Of course, this analysis relies strongly on the fact that 
no dynamical effect save the rotation is present in the prob- 
lem. The major corrections to take into account would be 
turbulent motion, magnetic field and chemical gradient. In 
regions where turbulent motions are expected (i.e. possibly 
near the edge of a convective zone, or in regions of strong 
shear) their effect could in principle be included in equation 
(1), provided a Reynolds-stress tensor for the microscopic 
motions is known. However, it is likely that the turbulent 
Prandtl number is much closer to unity, which would then 
render the analytical analysis obsolete (but not the numeri- 
cal abilities of the code). If magnetic fields are present, either 
on large or small lengthscales, the whole dynamical struc- 
ture of the star may be changed completely. A good example 
of this effect is well illustrated in the sun; observations reveal 

t although only if magnetic fields and turbulent motions are too 
weak to influence that equilibrium 



that the radiative interior is rotating nearly uniformly, as the 
shear imposed by the convection zone is suppressed within 
a very thin layer (the tachocline) . Models exist (e.g. Gough 
& Mclntyre, 1998) in which a small-amplitude, large-scale, 
fossil magnetic field within the radiative zone is able to im- 
pose uniform rotation throughout most of the interior. This 
type of fossil field is not unrealistic in any star. Finally, gra- 
dients in the mean molecular weight may also significantly 
alter the dynamics of the system by reducing or augmenting 
the effective buoyancy of the fluid, especially near bound- 
aries with convective regions. Mestel (1953) has described 
how this effect may counterbalance the slow, rotationally 
driven motions, and Theado & Vauclair (2001) have argued 
that this effect may be essential in explaining the observed 
Li abundances in stars. 

As a result, the rotation profile in rapidly rotating stars 
may not resemble that predicted in Section [| but any ob- 
served discrepancy with the predicted model will be a clear 
signature of the effects discussed above. 
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